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Electron-electron interactions in half-filled high Landau levels in two-dimensional electron gases 
in a strong perpendicular magnetic field can lead to states with anisotropic longitudinal resistance. 
, This longitudinal resitance is generally believed to arise from broken rotational invariance, which is 

^SJ ' indicated by charge density wave (CDW) order in Hartree-Fock calculations. We use the Hartree- 

Fock approximation to study the influence of externally tuned Landau level mixing on the formation 
of interaction induced states that break rotational invariance in two-dimensional electron and hole 
systems. We focus on the situation when there are two non-interacting states in the vicinity of the 
Fermi level and construct a Landau theory to study coupled charge density wave order that can 
occur as interactions are tuned and the filling or mixing are varied. We examine in detail a specific 
example where mixing is tuned externally through Rashba spin-orbit coupling. We calculate the 
phase diagram and find the possibility of ordering involving coupled striped or triangular charge 
D I density waves in the two levels. Our results may be relevant to recent transport experiments on 

I , quantum Hall nematics in which Landau-level mixing plays an important role. 
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PACS numbers: 73.43.-f, 73.20.Qt, 73.43.Nq 



I. INTRODUCTION 

C ' The integer^ and fractional^ quantum Hall effects have been studied extensively in two-dimensional electron systems 
^ (2DES).'^ It has been proposed^i^ and observed^ that the states at half odd integer filling fraction in low Landau levels 
I I (LLs) have distinct behaviour from states at odd denominator filling fractions. This distinction extends to higher 
LLs, where anisotropic transport has been observed at half odd integer filling fractions/"^ which is believed to arise 
' from broken rotational symmetry due to electron-electron interactions. ''^^^^'^ The experimental evidence for broken 
rotation symmetry at half odd integer filling fractions is large peaks in the longitudinal resistivity, which vanish as 
the direction of applied current is rotated by ninety degrees. 
', Hartree-Fock calculations of the effects of electron-electron interactions in high Landau levels indicated a transition 
' to a striped charge density wave (CDW) phase at half odd integer filling fractionsi^i"— It has been argued by Fradkin 
l' [ and Kivelson and others that the anisotropy is evidence of an electronic nematic phasei ^^d^" — However there continues 
■ to be debate as to whether the states are electronic nematics or smecticsi^"— 

' There have been several recent experiments on two-dimensional electron or hole systems that are in or close to 
the regime where anisotropic quantum Hall states are expected for which Landau level mixing is believed to be 
either appreciable or tunable via an external parameter i^""— LL mixing can arise from a variety of sources, such as 
disorder, geometric confinement effects in a quantum well, interactions or spin-orbit coupling^^ but has generally been 
' ignored in theoretical studies of anisotropic quantum Hall states jiSrJ^iii A couple of exceptions are investigations of 
H I the effects of LL mixing from electron-electron interactions^^ or disorder— on CDW formation which may be relevant 



for experiments showing re-entrant behaviour of quantum Hall states in the lowest LLi^ There has also been recent 
interest in LL mixing in the context of non-Abclian quantum Hall statesi^"— 

We investigate the effects of Landau level mixing on the formation of quantum Hall states with broken rotation 
symmetry by using the Hartree-Fock approximation to study CDW formation in a two-dimensional system of charged 
fermions (electrons or holes) when LL mixing and interactions are present. This neglects quantum and thermal 
fiuctuations,^^ but can be hoped to be as informative as Hartree-Fock studies on single Landau levels . '^^'^^ We assume 
that there is a term in the Hamiltonian that mixes LLs, leading to non-interacting single-particle eigenstates that are 
linear combinations of LLs. We explore the effects of interactions on CDW ordering within these eigenstates. We 
specialize to the situation when two energy levels lie close to the Fermi energy, Ep, of the system, which can lead 
to competition between CDW phases of various symmetries originating in different non-interacting states. To study 
this, we construct a Landau theory for CDW ordering in the presence of LL mixing when there are two states that are 
close to the Fermi energy. In the range of energies where the two states are close to degenerate we also consider the 
possibility of quantum Hall ferromagnetism, but for the example we consider here, CDW ordering occurs at higher 
temperatures within the Hartree-Fock approximation. 

In order to study the effects of LL mixing on CDW formation experimentally in a systematic way, some external 
tuning parameter that can mix LLs is required. Neither disorder nor interactions are particularly suitable for this 
task. However, Rashba spin-orbit coupling can be tuned experimentally in a quantum well by parameters such as 
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applied gate voltage. This motivates us to investigate the phase diagram as a function of filling fraction, LL mixing 
and temperature when LLs are mixed via Rashba spin-orbit coupling when two states are relatively close to the Fermi 
level. We find that at half-integer filling fractions there can be triangular CDW ordering in both levels simultaneously 
when the filling of the higher energy level is non-zero. 

This paper is organized as follows. In Sec.|lT]we investigate electron-electron scattering in the presence of LL mixing 
in high LLs, using a similar approach to Moessner and Chalkerjl^ to argue that there is still an instability towards 
CDW ordering in the presence of LL mixing. In Sec. IIIII we construct a Landau theory for striped and triangular 
CDW order when there are two states close to the Fermi energy. In Sec. IIVI we display numerically calculated phase 
diagrams for LL mixing due to Rashba spin orbit coupling and also discuss possible quantum Hall ferromagnetism 
near degeneracy. We discuss our results and their possible connection to experiment in Sec. [V] 



II. INTERACTIONS IN THE PRESENCE OF LANDAU LEVEL MIXING 

We consider fermions (electrons or holes) confined to two dimensions in a perpendicular magnetic field. We label 
single-particle basis states as \nkjz), where n labels LL index, k is the pseudomomentum, and jz is the spin eigenvalue. 
The projection of the orbital part of the eigenstates onto real-space is the LL wavefunction 



^i") (r) ^ _=4_e^'^-^0„ ( ] , (1) 



where we work in the Landau gauge A (r) ~ B {—y, 0, 0), with L the spatial extent of the system and Iq = y/h/eB 
the magnetic length. 

We assume that there is some source of mixing of the basis states, such as spin-orbit coupling, which leads to a set 
of non-interacting single-particle eigenstates of the form 

m = ^*"|afc), (2) 

a={"Q J?} 

where the C$a are the coefficients of the basis states and \ak) is shorthand for \nakj"). 

We now allow for interactions between fermions via a potential V (r), which has Fourier transform V (q), leading 
to the interaction Hamiltonian 



X 



j2„ 



2 ^ '^*ifci"*2fc2"*3fc3"*4fc4 ^ C'|.iqC'|.2CT^*3pC*4^'^j 

{*1,'I>2,'I'3,*4} {fcl,fc2,fc3,fe4} a/Sap 
,2 

F(q)(n„fci |e-'i--i|n^fc4)(n.fc2|e^'i--^|npfc3), (3) 



{2nf 



where a^k is the annihilation operator of \^k), and the four single particle states involved in the interaction are 
labelled {<I>i, $2, ^3, $4} with corresponding pseudomomenta {fci, ^2, k^, ki}. The orbital matrix elements have been 
evaluated by Raikh and Shahbazyan^ to be 

(4) 

where 

\ nali\ I "a/3 
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with nai3 — niax(na,n^), n^^ — min (rtQ,, n^), Arto,^ = \na — n/^l, and 9 — aictan (qy/qx)- Equation ^ can be 
rewritten as 



H 



{*1 ,*2,'E'3,'E'4} fc'™ 



'l>2,'l>3,'i'. 



where 



and 



-^{*i,*2,'I'3,*4} ("^1 — E ^*ia^*2<T^*3pC'*4/9^a^^crp-^acrp/3 (JTI, /) , 



(5) 



(6) 
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The function M^^^ $2.'i'3,*4} C*^: describes the strength of the interaction between the eigenstates $1, $2, 'i'a, and 
$4 (ihustrated in Fig. [1]), with the function laapp {mj) describing the overlap between the orbital part of the basis 
states. The strength of the interactions between particular eigenstates can be tuned by using external parameters to 
vary the coefficients, C$.q. 



$1 , k + I 



$2, k + I + m 





$4, k $3, k + m 

FIG. 1: Interaction between states <E>i, $2, $3, and $4. k, I, and m are pseudomomenta. 

A. Hartree-Fock Approximation 

We now investigate the self-energy and two particle scattering vertex using the Hartree-Fock approximation. Moess- 
ner and Chalker^^ showed that for fermions interacting with the hard-core potential V{q) — —vq^ within a single LL, 
the Hartree-Fock approximation is exact in the limit of large LL index. We generalize this result to electrons with 
different LL indices interacting via a hard-core interaction. Provided the difference in LL index is much less than 
the maximum LL index, then in the limit that the maximum LL index is large, the Hartree-Fock approximation is 
asymptotically exact for the hard-core potential. We return to consider Coulomb interactions in Sec. IHII 

To take the large-LL limit, we follow Moessner and Chalker and replace An^nf,{x) by its WKB envelope: 
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The envelope is defined for x"_^ < x < x"f , where 
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We now turn to the properties of the interaction vertex -^^{$i,$2,'i'3.*4} limit of large LLs; specifically, we 

study the properties of function laapp ("t-, [Eq- ©]■ The function laapp {'m, gives the strength of the interaction 
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between eigenstates in terms of the constituent basis states. For the single-particle and two-particle vertices, we focus 
on the situation when there are two eigenstates lying near the Fermi level. 
Substitution of Eq. ^ into Eq. (O yields 



laapll {m, I) = — 



("a/3 + 



1/2 



7,' I 



1/2 



dq, 
27r 



V {qx,l) cos /o9xm - {ua - + ncr - rip) (arctan - |^ 



, An„ 



1/4 



, An„ 



-,1/4- 



(9) 



As the eigenstates in the presence of LL mixing are not Landau levels, when we refer to the large-LL limit, we mean 
that the non-interacting single-particle eigenstates near the Fermi energy are dominated by contributions from Landau 
levels with large indices. If Auap ^ 0, and AtIq.^ <C n'^p, then 



n' I 



«/3 + Ana,3) 



a/3 



as n'^p oo, which causes the interaction vertex lacpp (m, I) to vanish. Thus, in this limit, laapp (jn, I) is finite only 
when AnQ« — ts.n„p = 0, in which case 



' Jo 



dqx 



V{qx,l)cos [loq^m] 



2tt 



v/^fTP [8n, - {ql + [8n. - (g^ + P)] 



72^1/4' 
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where 



max 



^{q-^-^f-P, 
max { \/%na, \/8rv } • 



If Ua = Ucr, the envelopes in Eq. (jlOp associated with Ua and overlap exactly, whereas when Ha ^ n„ there is less 
than perfect overlap, so laaaaim, I) will be an upper bound on the interaction vertex. 

Let N denote the index of the highest LL basis state that contributes to occupied non-interacting eigenstates. 
Moessner and Chalker showed that for a contact potential laaaaimjl) is non-zero only when m = in the large-LL 
limit, hence, as we have established that laaaaim, I) is an upper bound for Iaapi3 {m, I), it follows that Iacrpi3 (w, I) is 
nonzero only when m = 0, and obeys the bound 
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Thus, in the large-LL limit, only exchange terms are important in the interactions between basis states for a contact 
potential. The results of MC for other classes of diagrams can also be applied here as an upper bound to other types 
of diagrammatic contributions. Thus, diagrams with crossed interaction lines and closed fermion loops vanish in the 
large-LL limit. When we consider the case where the relevant constituent LL indices are not necessarily large, we 
apply the Hartree-Fock approximation directly. In this situation, both Hartree terms and Fock terms are important. 

To avoid divergence of the bound on laapp as the number of constituent basis states increases, we scale v with M 
such that vM is constant. Assuming that each eigenstate $ is a linear combination of P LL basis states (for example 
for Rashba spin-orbit coupling considered in Sec. IIVI P = 4), we obtain the following upper bound on the interaction 
vertex: 



2vVI\f 
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1. Smgle-particle propagator: Self-energy 

We now turn to consider the self-energy in the presence of LL mixing. The finite-temperature Green's function for 
an eigenstate $ is 



where a;„ = (2n + 1) tt/ [3 is the n"^ Matsubara frequency, = — fi, where is the energy of non-interacting 
eigenstate $, /i is the chemical potential, and E,i> is the self-energy. We note that the bare Green's function, 
G'o($,a;„) = (iw„ — is independent of pseudomomentum k: the LL mixing is chosen so that the energy is 

independent of fc, as for single LLs. 

In this section we determine I]$ using the Dyson equation, keeping only terms which do not vanish in the large-LL 
limit. Note that all Hartree diagrams vanish for the special choice of a contact potential, so we need only consider the 
Fock diagrams. For a diagram with a single propagator and a single interaction line, the summations over Matsubara 
frequencies and momenta / decouple. Then 

- ^ G($„,w^)M$$„$$„ (0,j). (11) 

The interaction vertex corresponding to (Oij) is illustrated in Fig. [TJ Since we consider only exchange 

interactions in Eq. ()11|) . we adopt the notation Af$$^$(i)^ (0,i) = Af$$^ (Oij) for brevity. 
The full Green's function can be expressed in terms of the bare propagator. Go ($, a;„) as 



G($,w„) = Go($,cj„) + ^Go($,w„)'^ ^ Go($„,^m)A^$*„ (0,j) 

+ ^Go($,w„)' Go($c«,^m)M3,$„(0,i)Go($^,Wp)Af*„$, (0,fc) 



+ ... . 

We can simplify each term in the expansion by noting that the sums over Matsubara frequencies and momenta 
decouple. Let A^$^$^ — (Oii)- Then 

S$ = -^^^Go ($a,W™) A^cE,*^ - ^ ^ Go ($a,Wn)^ A^cE,4,„Go ($^,W™)7W$„<1>^ 

X! X! Gq ($a,W„)^ A^<E.<E,„Go ('I>/3,W„0^ A^<I.„<E.^Go ($^,Wp) A^<t,^<[,^. 

<E>Q ,<E>^ ,<I>^ ii.Jmi'^p 

Noting i J:^^ Go (<f>,^^„)=/ ($$ - m), ^ Ec.„ Go (<f>, t^„)'=-/3/ (C$ - m) [1 - / (e$ - m)], and (e$ - m), we get 

^ (1 - i^*„) A^*<i.„i^$3 (l - i'*^) A^*„*^i^*^A^$^<i.., + ■ • ■ • (12) 

,<E>^ j*.^ 

^. Two-particle vertex: CDW InstabtUty 

We now turn to investigate the two-particle vertex, which was shown by Moessner and Chalker to indicate an 
instability towards charge density wave (CDW) order in the case of a single Landau level. We show that such an 
instability persists when two levels are involved and identify temperature scales corresponding to CDW ordering in 
the specific example when there are only two states, $i and $2: near the Fermi energy. Moreover, we identify a third 
temperature scale related to the possibility of mixed order, which we discuss in more detail in Sec. [ 



6 



*i , k + q 



*4, k 

w + A 




$3, k + q 

n 



4)2, k 
A 




FIG. 2: Diagrammatic representation of Bethe-Salpeter equation. 



Consider the two-particle vertex, which can be determined via the Bethe-Salpeter equation (Fig. [5]): 

r$„$^$.^$j — Af$^$^$^$_5 (0, g) 

(13) 

The summations run over all momenta, Matsubara frequencies, and eigenstates. 

First, we examine the case when uj =^ 0. We rewrite the sum over states as two separate sums: a single sum over 
for the terms where $j ~ along with a sum over $j ^ The first sum vanishes when u ^ 0, whereas for the 

second sum, —i J2x {^k,^ + A) Go (<&j, A) — - — jj — ^ — r. Then, the two-particle vertex is 



(<Z, — (0, 

2^A^*fc*f,*j*i (0, Or*„*j*^$fc [q- 7- 



(14) 



Now, consider uj — We separate the summation over eigenstates into two cases: "I>j = $fc and $j 7^ If 
- then i Go A)' = -^v^^ (l - . If ^ then i Q ($,, A) Gq ($fc. A) = ™" 
leads to 

r$„$^$^$^ (g, 0) = M$^$^$^$j (0, g) + (l - 1^*3) Af$^.$^$^$^ (0, Z) r$^$^.$^$^. (g — Z, 0) 

_ ^ V^ii — ^M<i>,,<[.^$^.$, (0,Or<i,„<i.^$^$, (g-/,0). (15) 

Equations (|14p and psp are exact expressions for the two-particle vertex at finite and zero frequency, respectively. 

The expressions in Eqs. (|14p and (jlSp are somewhat cumbersome. Hence, it is instructive to investigate a simpler 
case than the general situation, in particular, the case where only two states are located near the Fermi energy. 



3. Two states near the Fermi energy 



Let $1 and $2, with energies < ^$3, label the two states that lie nearest the Fermi surface and assume that 
both and z^^^ are non-zero. We will ignore other states that are either fully occupied or empty. Furthermore, 
orthogonality of $1 and $2 leads to the simplification that only four vertices contribute to the one- and two-particle 
propagators: M$^$j$j$j (0,g), M^j^^^^^^ (0, g), M^^^^^^^^ (0, g), and Af$2$i$i$2 (0,g). 

The form of the vertices that are non-vanishing implies that A^*„$,j = ■^$Q*a<5$^$^ , so from Eq. (fT2|) the self- 
energy of the state $ is 

* 1 — jSv^ (1 — i/cj)) A^$* 
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Since Ad is negative, the self-energy has no poles, and hence is well-behaved for a contact potential. 

We can also obtain simple expressions for the two- particle propagator. The condition of orthogonality of $i and 
$2 limits the possible non-zero terms to r$j$j$j$j, r$j^$2$2$j, T^^^-^^-^^^, and r^^^j^^^^. First, consider the case 
of finite frequency. Taking the Fourier transform of each expression with respect to q with conjugate variable s, we 
can solve for each of the two-particle propagators explicitly. We find that r$^$^$^$^ = M$j$j$j$j (0, s), and 

r$24'2*2*2 i^j^) — Af$2*2*2*2 (Oj*) wfilch havc no poles, and 

i<i,i<i,2<i>2<i>i (s,a;) - —J- v^,-u:^^ ^ (l^J 

1 + A/*i*2*2*l (0, s) -—^^—^ 

^*2<E'i*i*2 - 77— ^*2-''*i ■ ^^'> 

i - Ai$2$i$i$2 (L), S) —7- r 

The r$j$2*2*i ('^1 and r$2(i.j(i.j(i.2 (s, w) propagators have poles at finite w when 7^ j/$2 ■ Since there are no poles 
in the Re (w) > half-plane, we can analytically continue Matsubara frequencies to real frequencies, io; — > -f i5. By 
symmetry, M$j$2*2*i (Oi ■s) = -^$2*1*1*2 i^^^) Eqs. and ([T7)) indicate that the system is most susceptible to 
ordering in states $1 and $2 at the same wavelength scdw 

Let Az/ = i/$2 ~ and = ^$2 ~ : and assume ^$2 > C*i • After analytical continuation, the frequency at 
which r diverges is 



oJc = A^ 



1 + M$^$2*2*i (0; 



At low temperatures, uJc oc /3, which implies that uj diverges as T — > unless w — > at finite temperature. This 
motivates us to consider the static {uj — 0) limit. 

In the static limit, all four two-particle vertices are modified by the interaction. The vertices are 

f, / a\ Af$i$i$i$i (0,s) 

1*1*1*1*1(3,0) = 7——, (18) 

/ n\ Af$2*2*2*2 (Oj /1n^ 
1*2*2*2*2(5,0) - -; T^r-T' (.19) 

1 - Pm-i (1 - 1^*2) ^■^*2*2*2*2 (0, S) 

and 

f, /' — Af*i*2*2*i (0, s) /on^ 

1*1*2*2*1 (s,u) - -, , ^^2-^-1 A> ?;r~:' ^ ^ 

1 + cJ7-^^^*i*2*2*i (0> s) 

and by symmetry of the interaction vertices, r$j$2*2*i (■^lO) = r*2*i*i*2 ■ We see in Eqs. and (IT^ that 
the system is unstable to CDW formation in $1 at a temperature 

Til — >^<s>i (1 — i^*;) A/$i<i.i<i.i<i.i (0, si) 

and in $2 at a temperature 

^22 — ViS>2 (1 — ^^$2) Af$2*2*2*2 (0i'52) , 

where Si and S2 maximize Af$i$i$i$i (0, s) and Af$2*2*2*2 (Oj respectively. 

If f $2 ~ '?*i is small, then the filling fractions j^$i and are almost equal. Expanding v^^ about to first order, 
1^*2 ~ V'S'i + /' (^*i - m) (C*2 - ^*i)- Then ^£2^^ « (1 - z/*, ) + O (/"' (Ci)) • Eq. (I20l) allows us to find a 

third temperature scale, associated with interactions between states 1 and 2: 

ri2 = 1^*1 (1 — Z^*i) Af$j$2*2*l (0i'5l2) , 

where S12 maximizes Af$j$2*2*i (O7 Whilst we have derived these results for a contact potential, we will see that 
multiple temperature scales also arise for a Coulomb potential in Sec. IIVI 
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III. LANDAU THEORY FOR MIXED CDW ORDERING 



There has been considerable work on the Landau theory for interaction induced CDWs in the lowest Landau 
level,— 1^ and higher Landau levels j ^ '- l ^ Here we allow for ordering of CDWs of states that arise from Landau level 
mixing and focus on the situation in which there are two states close to the Fermi energy. We use similar notation to 
that of Fukuyama et al. (FPA)J^ Furthermore, we modify the notation of previous sections and write J2<s> ~^ Si- 
The analysis in this section leads to a Landau theory where coupling between order parameters of different states 
becomes important. 



A. Hartree-Fock Hamiltonian 



We consider an interaction u(q) and start from the interaction Hamiltonian 



(21) 



which is a natural generalization of the FPA interaction Hamiltonian to a situation with multiple statesi^ The 
operator p (q) is the charge density operator of the two-state system, expressed as 



/5(q) = J2J2P'^ ^'^''"'Ix+'^^x. 



i=l X 



where 



P^x{cl) = y'd2re-"i-V*x(r)<^«(r)=e-^''-^A,,(q), 
where (/)ix(r) = J2a C'iaV'x °^ with defined in Eq. ([1]), and with 



where n'^ 
operator 



-«("ic.-"ja)(^-'^/2)f 



Sriij /2 



(qlo) 



(22) 



(23) 



(24) 



I, and 9 = arctan (qy/qx). The operator Sx (q) is a pseudospin density 

(25) 



-S':.(q) = J^Jl^'i'^^^^^-iX+ajx 

= l X 



where 



%,x(q) = al^ Jd're-^^-'^(t>*x(^^^x{r) = al^e-^'>^''A,{q), 



(26) 



and ax is the Pauli matrix. The first term in Hj describes interactions between fermions in the same eigenstate, 
while the second term describes interactions between fermions in differing eigenstates. Such a term can generically be 
expected to be present in the situation we consider herej-3*. but usually will have a smaller magnitude than the first 
term in Hj (which is the case here). The final term in Hj is 



/(q) = Y.Aui<lfpm- 



(27) 



We use different notation from the notation in Sec.|TT]by labeling states with X = kl^ rather than pseudomomentum fc, 
and X± — X± ^qyl^- The indices i and j take values of 1 or 2 and label states $i and $j. The operators a-j^-^ (aix± ) 

create (annihilate) a fermion in the non- interacting eigenstate $i with guiding centre coordinate X± = X ± , and 
satisfy the usual anticommutation relations | <iix±, '^jY± f ^ ^ij^x±,Y±- In the limit of no mixing, and when restricted 
to the lowest LL, the interaction Hamiltonian reduces to the form considered by FPA. 
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The total number of fermions in the two-state system is N = Ni+N2^ where 



(28) 



X 



is the number of electrons in state $i. Similarly to FPA, we define the order parameter for a CDW in eigenstate 
A, (Q,) by 



Qi 



Qy 



(29) 



where Pi = Ni/N is the fraction of the fermions that occupy state $i. Each order parameter depends on a unique 

set of wavevectors {Qi}j where Qi = (yiQi)^ , {Qi)^- The summation over wavevectors allows for CDW ordering at 

several wavevectors. For instance, a triangular CDW can be described by a set of three wavevectors of equal length, 
oriented at angles ^ with respect to each otherfiS, whereas a striped CDW is described by a single wavevectorJ^ 
Furthermore, note that the order parameters are defined such that A* (Qi) ~ A^ (— Qi), which follows from Eq. p9p . 
Next, we use the Hartree-Fock approximation to rewrite the Hamiltonian as: 



where 



H = Ha + Hcuw + ^E E (^') '^(^ 

«=i {Q.} 



Ho = ^^<^^alxaix, 

i=l X 



is the non-interacting Hamiltonian with the energy of single particle eigenstate <J>i, and 



(30) 



HcBW = - < E E (^'') E (Qi) e'^'^'^''''4x.a^x^ + A* (Qi) e-^(«')-^al^^a,x_ 

I 1=1 [Q.} X 



-EE^^^'^-(^?^)E h (Q0e'^'5'^-^a]x_a.x+ +A:(QOe- 



"■jX+"-jX- 



i=l j^i 



(31) 



We let Qi = IQil and 9i ~ arctan and henceforth we absorb Zq in the definition of q, i.e. g/g ~^ 1- In the 

summations over wavevectors, we sum over {Qfj^^^' '^h^re n is the number of wavevectors describing the CDW 
state. Furthermore, we define the following functions, which describe the Hartree-Fock interaction potentials: 



Uii {Qi) = I dqqv{q) Ai^ {qf Jq {qQi) - v {Qi) Ai {Qi)^ , 



1 °° /■ 

U^]{Q^) = ^ E N^c^v{q)A,,{q)Aj^{q)J„,{qQi)~v{Q{)A,{Qi)Aj{Q,), 

m— — C30 

The potentials J7ii and Uij represent the most significant contributions to the Hartree-Fock Hamiltonian. We ignore 
other terms generated by the Hartree-Fock approximation, as their contributions are negligible. 



B. Landau theory 



We construct a Landau theory in the spirit of FPA by considering the difference in free energy between a uniform 
system and an ordered system for a fixed total number of electrons, N . We write the difference in free energy, 5F ^ in 
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terms of the grand potentials of each system, 

where rJcDW a-nd f2o a-r^ the grand potentials of the CDW and uniform systems, respectively, and the chemical 
potentials of the uniform and CDW systems are fiQ and /i, respectively, and we expand in powers of the CDW order 
parameters Ai (Qi) and A2 (Q2)- 

The grand potential when there is CDW ordering is 



(fi) = -/3-iln|Tr [e-^(^-^^-^^"'^)]} 



CDW 

= no (/i) - r ' In (e-^^^°»)„ + E ^'f^- (Q^) (Q^)l' ' 

i=l {Q.} 

where (. . .)q indicates an average with respect to the uniform state, Hqdw is defined above in Eq. pip , and 

i 

We now develop the perturbative expansion of the free energy in terms of the order parameters Ai (Qi) and 
A2 (Q2)- We expand the logarithm, which gives 

The p = 1 term in the expansion vanishes, since it involves averages of the form (al-^_^aix_ ^ = Sx^X- which imply 

Qi = 0. We consider Eq. ([32]) to fourth order in both order parameters and write the CDW order parameters in the 
form Afe = Afc (Qfc). The general form of the free energy density when Qi 7^ Q2 is 

:f ^ Plai (Qi,r) |Ai|2 + P|a2 (Q2,T) + PfW {Qi,T) |Ai|3 + P|62 (Qi,r) |A2|3 

+^1^01 (Qi,r) |Ai|4 + Pic2 (Q2,T) |A2|4 + PfP|7(Qi,Q2,^,r) |Ai|2|A2p, (33) 

where J- ~ 5F/N. We calculate the coefficients a^, 6^, Ci and 7 explicitly within the Hartree-Fock approximation 
allowing for the possibility of either striped or triangular CDW order in either state 1 or state 2. The striped phase 
is described by a single wavevector while the triangular CDW is described by a set of three wavevectors of equal 
magnitude, oriented at 120° with respect to each other. The expressions for these coefficients are quite lengthy and 
we present them in full detail in Appendix [X] The phase behaviour of free energies of the form in Eq. (1331) when 
^1 = ^2 = was studied in detail by Imry.'^^ In general, there are three possible critical points, with free energies 
J^i, J^2j and Tm- If is globally minimum, then Ai ^ 0, and A2 — and vice versa for J-2. If Tm is the global 
minimum, then both Ai ^ and A2 ^ 0. In the case we consider, the phases associated with Fi and F2 are CDW 
ground states in either state or $2, respectively. The phases associated with Fm are mixed CDW states. The 
mixed CDW phase exists only when the parameter 7, which couples ordering in states 1 and 2 is finite. If 7 > 0, then 
there is competition between order parameters of each state; if 7 < 0, then ordering in one state serves to enhance 
ordering in the other The phases of the order parameters are important only for triangular CDWs If the phases 
corresponding to the three different order parameters are denoted by {ipiYl^i, then for a state $, + 02 + 03 = (tt) 
if i^<E> < 1/2 (> 1/2). 

There are additional terms in the free energy that are allowed when Qi = Q2, which take the form 

^•^Qi=Q2 = PiP2a|Ai||A2|+Pi2p2?7|Ai|'|A2|+PiP2V|Ai||A2|'+Pi'P|7|Ai|'|A2|' 

+PfP20 |Ai|3 IA2I + PiP|p |Ai| IA2I' . (34) 

The coefficients a, rj, Lp, 7, ct, and p depend on the magnitude of the ordering wavevectors Qi — Q2 — Q and the 
relative phases of the order parameters. Note that 7 modifies 7(Q,Q,0,T) in the case of equal wavevectors. These 
coefficients are displayed in Appendix [Bl 

We adopt the following notation for describing the possible phases. For ordering in a single state, we label according 
to the symmetry and state index. For example, striped (triangular) ordering in $1 is denoted SI (Tl). For mixed 
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CDW ordering, we label according to the symmetry of the order parameter for each phase with order in $i preceding 
that in $2- For example, mixed CDW ordering between striped phases in both states is labeled SS. If the wavevectors 
for each phase are equal, we add an extra "E" to the end of the label. In all, we allow for the possibility of 12 types of 
CDW ordering: SI, S2, Tl, T2, SS, TT, ST, TS, SSE, TTE, STE, and TSE. We present results for specific numerical 
examples in Sec. IIVI 

IV. NUMERICS 

The Landau theory developed in Sec. IIIII makes no assumptions about the source of LL mixing, and does not specify 
the interaction v(q). There are many possible sources of LL mixing in 2DEGs that can affect CDW order, such as 
electron-electron interactions'^^ or disorder. "^^ However, we focus on the specific example of Rashba spin-orbit coupling, 
since this may give a way to tune levels into degeneracy controllably and create a situation of the form discussed in 
Sec. IIIII Electron-electron interactions are described by an unscreened Coulomb potential ~ we expect that screening 
may lead to small quantitative changes in our results, but that our qualitative results will be robust as we discuss in 
Sec. lIVA2l 

A. Model and parameters 

1. Rashba spin-orbit coupling 

Rashba spin orbit coupling can be significant in GaAs quantum wells and leads to LL mixing. For holes near the F 
point in GaAs there are four j — 3/2 spin states in the valence band, which separate into two heavy- hole (HH) and 
two light-hole (LH) states with different effective masses. The heavy hole states correspond to angular momentum 
jz — ±3/2 and the light hole states correspond to jz = ±1/2. The strength of Rashba coupling can be tuned through 
external means, such as a back gate, which allows tuning of levels into and out of degeneracy. Generically, Dresselhaus 
spin-orbit coupling will also be present in GaAs, however, we focus on the regime where Rashba coupling is dominant 
and ignore additional (non-tunable) mixing that may arise from Dresselhaus coupling. 

We specify the single-particle Hamiltonian, which includes Rashba coupling and Zeeman splitting, as 

-ffnon-intcracting = Hq+Hz+HsO, (35) 

and diagonalize to find the non-interacting eigenstates as a function of the strength of the magnetic field and the 
Rashba coupling. We then consider a specific example in which there are two energy levels that are close to degeneracy 
at a particular value of Rashba coupling, and explore the effects of interactions on the phase diagram as a function of 
temperature, filling and Rashba coupling. 
We take the free-particle Hamiltonian to be 

Ho = Y,fu^{\jz\) (fij^+l 

where 



±f 
±^ ' 

and the number operator fij^ = a^^ , with the LL raising operator for a hole with Jz eigenvalue jz ■ ojhh and 
LJLH are the cyclotron frequencies of the heavy and light hole states, respectively. The Zeeman splitting term is given 
by 

Hz = -g*fJ.BBJz, 

where the effective g-factor for the holes is g* . The third operator may be written as 

Hso = -a7^8^,8^, °) (k X J)^ , 

which introduces Rashba spin-orbit coupling into the Hamiltonian. The parameter TZ^vSv is a material-dependent 
Rashba coefficient; {£z} — is the electric field due to the fermions in the absence of a substrate fieldi^ We use 



UJHH for jz 
^^LH for jz = 
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parameters that are closely related to those relevant to Fischer et al.'s experiments. We choose material parameters 
g* = 1.2, TZsvSy — 14.62 eA , and e = 12.4 eq, appropriate to bulk GaAsj^ we set the heavy- and light-hole masses 
equal to tojj = 0.51mo and to^ = O.OSmo, respectively)^ and we set the in-plane carrier density to p = 3.5 x 10^"'^ 

— 9 41 

cm ^iSi 

We model the effect of tuning via external means (such as a back gate) by the parameter a. The magnetic field 
B is fixed according to the total filling fraction of the system, vt = l-Kpl'^. The term a (f^) may be compared to 
the average confinement electric field seen by fermions in a triangular potential, as discussed by Fischer et al^ The 
parameter a then corresponds to the enhancement of the average electric field by increasing the substrate electric 
field, £s- That is, a = 1 -I- -^j&y- For no substrate field, a = 1, whereas in the Fischer et at. experiments, a ~ 2-3 
corresponding to nonzero 8s- For the levels we study here, a crossing takes place for aa ~ 20. The values of a at which 
level crossings occur are very sensitive to the model parameters. In particular, the values of the material parameters, 
such as the effective masses and the carrier density, and the presence or absence of Dresselhaus coupling can alter 
the positions of crossings and the nature of the levels that cross dramatically. Fischer et al. construct a Hamiltonian 
from the Fg, Fg, and F^ bands, whereas we consider the Fg band only. Second, they include both the Rashba and 
Dresselhaus spin-orbit coupling terms, whereas we consider Rashba coupling only. Third, we set the carrier density 
to the upper limit studied by Schmult et al.,'^^ which is larger than any of the densities studied by Grayson et al. 
This choice of carrier density increases the energy gap between the energy levels comprising the two-state system and 
the nearby energy levels, as shown in Fig. (jS]), thereby allowing us to focus on states $i and $2 only. Whilst our 
numerical calculations are not adequate for a quantitative discussion of the experiments in Ref. |2^, the crossings we 
are able to induce with increasing a in our simplified model should lead to a qualitative understanding of the effect 
of LL mixing on CDW formation. As we point out in Sec. IIV C 31 there are certain features in the phase diagram we 
expect to be generic and others that may depend more sensitively on details of LL mixing. 

We diagonalize the Hamiltonian shown in Eq. ([55)) and find the eigenvectors and associated eigenvalues as functions 
of a and magnetic field. The non-linear dependence of the energy levels on magnetic field allows for energy crossingsi^ 
Cubic and tetrahedral corrections can modify these crossings to anti-crossings.?^ Motivated by the results of Fischer 
et al^^ we focus on parameter values at which there are two energies that are close to degenerate and located near 
the Fermi energy. This is the situation we discussed in Sec. IIIIl In the vicinity of such a crossing, we write the 
filling fraction v = D + vq, where D is an integer, corresponding to the number of fully occupied energy levels, and 
I'o — 1^1 + V2 is the total filling fraction of the two partially filled energy levels. We label the crossing states $1 and 
$2, their associated energies Ei and i?2, and their filling fractions vi and V2. Our analysis is not greatly affected 
if crossings are replaced by anti-crossings, as the key point is that there be non-zero filling in the two states, which 
will be the case when the temperature is not too low for an anti-crossing. As T — > 0, only the lowest energy state is 
occupied, which decouples the two states, except when Ei = -£2- 



2. Interactions 



We choose the interactions between holes to be given by the unscreened Coulomb interaction, similarly to FPAiS 

27r e2 



v{q) = 



A 47re |q| ' 



where A is the area of the 2DHS (the factor 27r/yl comes from the Fourier transform of the real-space Coulomb 
potential). This allows for a somewhat more realistic interaction than the contact potential considered in Sec. HIl 

When there are multiple filled Landau levels, the filled levels can act as a background that screens the bare Coulomb 
interaction. It was found by Aleiner and Glazmar4^ (AG) that the screening from filled Landau levels can be accounted 
for by replacing e with 



2 



tAG{q) = <!l + — l-JoiqRL? S>eo, (36) 
I gas L J J 

where Rl = vp/ijJc and as are the Larmor radius and Bohr radius, respectively. Fogler, Koulakov and Shklovskiii^ 
used this screened interaction to show that electrons in a partially filled LL near half-filling preferentially order in 
a striped CDW phase in a weak magnetic field. Away from half-filling, they found ordering in a triangular CDW 
phase 

In the rt*^ LL at wavevectors gZo ~ I/v^j ^ag ~ RlIo-b\ at low magnetic fields, e^Q is large and screening effects 
are appreciable. However, at wavevectors gZo ^ 1 or ql^ < 1/f^, screening is negligible. In the example considered in 
this paper, the CDW ordering generally occurs when ql^ > 1 for an unscreened Coulomb interaction. This is outside 
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FIG. 3: Plot of energy of non-interacting eigenstates $1 and $2 versus tuning parameter a a,t ut = 7.5. Energy is in units 
of Eo = heB/nio and a is dimensionless. The nearest energy levels above and below Ei and E2 are labeled E4, and E3, 
respectively. Inset: Energy levels Ei and E2 for a £ [ao — 0.01, ao + 0.01]. 

the region where the effects of screening would be seen clearly if we considered the modified dielectric constant in 
Eq. ([36]). For the purpose of this work, we consider the bare Coulomb interaction as we focus on non- interacting 
states for which the dominant constituent LLs have n < 4, hence we expect our phase diagrams to be qualitatively 
correct, although screening may lead to small quantitative corrections. 

3. States 

Our approach to finding phase diagrams numerically starts by identifying a pair of states which have a crossing. 
In Fig. [3l we plot the energy of states that we label $1 and $2 for a magnetic field B corresponding to i^t = 8 
for our chosen electron density. In addition to these states, we also plot the nearest energy levels above and below 
$1 and ^2- At _B w 1.8 T the minimum energy gap to the levels above and below $1 and $2 is approximately 
0.23 meV ^ 2.7 K, which is large compared to the spacing between $1 and $2, which is approximately 2.3 /ieV. In 
comparison, the transition temperature we calculate for the CDW phase is ^ 10 K. The Hartree-Fock approximation 
thus overestimates the transition temperatures observed in experiment, in which anisotropic transport develops for 
T < 100 mKfiS, Nevertheless, we ignore adjacent energy levels and construct a two-state system from the crossing 
energy levels, and make only qualitative predictions about the possibility of CDW ordering as a function of vt and 
a. In the situation shown in Fig. |3] there are seven filled levels, and we investigate CDW ordering in the two states 
near the Fermi energy for 0.5 < I'd < 1. 

B. Quantum Hall Ferromagnetism 

We have focused primarily on the possibility of CDW ordering when two energy levels are in the vicinity of the Fermi 
energy, and each level has nonzero filling. When two levels are degenerate or close to degenerate then the physics 
of quantum Hall pseudospin ferromagnetism can also be important.— Following Jungwirth and Macdonald,'*'* we 
determine the Landau theory for quantum Hall pseudospin ferromagnetism when there is Landau level mixing using the 
Hartree-Fock approximation. This allows us to compare the ordering temperature for quantum Hall ferromagnetism 
with the temperature we have determined for CDW formation in the vicinity of degeneracy for the non-interacting 
states. 

We define the psuedospin operators in terms of the creation and annihilation operators of the non-interacting 
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eigenstates: 



(37) 



where aj is the Pauh spin matrix and I = 0, x, y, and z are pseudospin labels (cto is the two-dimensional identity 



operator). By construction, mj 



9y 



TO/., 



We focus on q = order, so the order parameter is defined as 



2n 



^" <; / t (y) 



(38) 



Since the total filling fraction of the two-state system is fixed, toq = i^t- 

We rewrite the creation and annihilation operators of states $i and $2 in terms of the pseudospin operators and 
proceed using the Hartree-Fock approximation in the same manner as outlined in Sec. IIIII or Ref. |4J. This leads to 
the effective Hamiltonian 



QHF 



(39) 



where 60 



1=0, z I,J=0,x,y,z 



IJ—0.x,y,z 



— \ Wii and hz — -^^j-^ — Wqz. The potentials Wij are defined explicitly in Appendix [Cl We 
note that for the eigenstates relevant to the numerical example shown below. Wo/ — Wzi — 0, I = x, y, and so toq 
and rriz do not couple to and my. 

To second order in the pseudospin order parameters, the change in free energy associated with ordering in a QHF 
takes the form: 



SfqHF = ^ bimi + ^ Wijmi (toj,o)o 

1=0, z I,J=0,z 



(40) 



IJ IK,JL 

The (••■)□ notation indicates that we average over non-interacting states. Thus, (too,o)o = vi + 1^2 ~ and {rhz.o)^ 
vi — 1^2 ~ Siy. It follows that 

r2 



SJqhf 



1, 



-WooTO^ 



! + (/{+ /2) (woo + + {f[ - /^) W,z 



+ lWyyml [1 - /? (l^T - 21^11^2) Wyy] 



+ l^Wzzml 



! + (/( + f'l) Wzz + 



Wzz 



+ {f'l - /2) W^z 



+ [vtWoo + SvWoz - bo] Too + [vtWqz + SvWzz + bz] mz 

+ [Woz {f[ + /^) (W^oo + Uzz) + (/{ - /^) (W^ooC/.. + Wl)] m^mz. 



(41) 



Since too is constant, dSfgHp/dmo — 0. Thus, terms containing too only act to shift the minimum of the free energy, 
and can be ignored safely. Let us then write the remaining terms as follows: 



SfgHF 



- E 

I—x,y,z 



ai (T) Wiimj + bz (T) rUz 



The coefficients a/ determine the temperature scales for a second-order phase transition into a quantum Hall ferro- 
magnetic state, whilst the coefficient bz acts as a symmetry-breaking field. We include the transition temperatures 
Tqhf obtained from the condition a/ = in the numerical examples in Sec. II V Cl below. 



C. Phase Diagrams 



We study the phase diagram as a function of a, ut and T in two different directions in parameter space. First we 
vary vt from 7.5 to 8.0 when the energy levels are degenerate. Second, we vary a at fixed vt = 7.5 and 8.0. It should 
be notedi^ that the free energy expansion is strictly only valid when the order parameters are small and is hence most 
applicable near the transition temperature into an ordered state. Thus, at a given a and vt, we consider the ground 
state to be the state with the highest ordering temperature. 
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1. Degenerate energy levels 

In Fig. 13] we show the phase diagram as a function of vt for the energy levels found in Sec. IIV A 31 when the two 
states are degenerate. As shown in the inset, when vt — 8.0, there is a transition from the uniform phase to the 
SSE phase. As vt is decreased from 8.0, the transition from uniform to CDW ordering prefers the TTE phase. 
Below the transition temperature, the transition between the SSE and TTE phases is first-order, and is a function of 
temperature and total filling fraction. Note that the Quantum Hall ferromagnetism temperature Tqh p is lower than 
the CDW ordering temperature for all ut considered. Thermal and quantum fiuctuations might change this ordering 
of transition temperatures. 

For each of the mixed states we find numerically that the order parameters Ai and A2 do not have the same 
phase. This can be understood by considering the lowest order term in the free energy that contains Ai and A2 (see 
Appendix [B|) : 

^•^Qi=Q2 = P1P2U12 iQ) [1 + Uii (g) /' {El ^ti) + U22 iQ) f {E2 - m)] IA1IIA2I cos (5012) . (42) 

Here, (5(/)i2 is the relative phase between the order parameters. The term in Eq. changes sign aX T — T12 — 
vi (1 — 1^1) Uii (Q12) + i^2 (1 — ^2) U22 {Qi2)- Since U12 {Q12) is negative for the values of Q12 that minimize the free 
energy near the transition, the sign of this term below T12 depends on the sign of cos((50i2). In the SSE state, 
5(j)i2 = TT, while for the TTE state, 02j - (/-h = tt, i = 1, 2, 3. Thus, the SSE and TTE CDW states can be thought 
of as coexisting CDW phases in $1 and $2 ordered at the same wavevector, but tt out of phase with one another. 

It is possible to determine another characteristic temperature associated with equal wavevector ordering in addition 
to T12. Truncating the free energy to O (A^) and minimizing with respect to both of the order parameters and their 
relative phases leads to the following quadratic equation whose solution yields a characteristic temperature Tm{Q): 

(1 - /3) r2 - 2 (ti2 - T + {T^^ - PT1T2) = 0, (43) 

where Ti = Un (Q) i^i (1 - i^i) + U12 {Qf jUxi (Q) V2 (1 - V2) (similar for T2), T^2 is defined from Eq. dH]), and 
/3 = 4[/ii (Q) U22 iQ) /U12 (Q)^- Solving Eq. (|43|) for T = Tm (Q) gives us an estimate of the transition temperature 
at which equal wavevector CDW ordering becomes relevant. Furthermore, maximizing (Q) with respect to Q gives 
us the critical value of the wavevector at the transition temperature. When terms beyond Eq. (|42)) are included in the 
free energy, it is not possible to obtain an analytical solution as we have above, and a numerical solution is required. 
T12 as defined above is larger than the ordering temperatures we find numerically and the presence of cubic terms in 
the free energy generically leads to a first order transition to the TTE state. 
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For the particular example we consider here, at the ordering wavevector Q, Uii{Q) /U22{Q) ^ 2, and we find that 
the ordering temperature for TTE is marginally higher than the ordering temperature for Tl which is much larger 
than the ordering temperature for T2. The tendency towards ordering is much stronger in state 1 than state 2, but 
the small energy gain from out of phase equal wavevector ordering in the two states implies that that the TTE phase 
has the highest ordering temperature. 



2. Phase diagram at fixed vt 



An orthogonal trajectory in phase space to that considered in Sec. IIV C H is to fix vt and then vary a and T. As a 
is varied, the filling fractions of the two states vary, which allows for phase transitions between CDW states at fixed 
temperature. Away from degeneracy, temperature also tunes the relative filling of the two states, with the filling of 
the higher energy state going to zero at low temperature. We plot phase diagrams for vt — 7.50 and 8.00. These are 
shown in Figs. [5] and |6] respectively. 
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FIG. 5: Phase diagram at vt ~ 7.5 in the a — T plane. Temperature is in units of ^^Jj^^,^ 
degeneracy point is indicated by a dotted black line on the a axis. 



and a is dimensionless. The 



At Vt = 7.50, we see that the TTE phase dominates around the degeneracy point. At vt — 8.00, the SSE state is 
favoured at values of a close to the degeneracy point, and the TTE phase is favoured away from degeneracy and near 
the transition temperature. Transitions between triangular and striped phases are first-order in nature. Away from 
degeneracy and at low enough temperatures, we expect only the state with lowest energy to have non-zero filling, 
leading to an integer quantum Hall state. The quantum Hall ferromagnetic ordering temperature is below the CDW 
ordering temperatures in both cases. Similarly to the degenerate case, ordering in state 1 appears to be the main 
driver of the transition with ordering in state 2 contributing a small lowering of the free energy. The balance of the 
contribution to ordering of the two states is non-generic and will depend on the details of the two non-interacting 
states that are participating in the ordering. 



3. Generic features of the phase diagram 



In our numerical example, we have focused on a particular pair of levels crossing at a particular filling fraction, so it 
is natural to ask to what extent the features we observe in the phase diagram are generic to when there is Landau level 
mixing as opposed to specific to the two individual levels that we selected. In general, the choice of two particular 
levels will affect the coefficients in the Landau theory through the mixing coefficients C$q (as introduced in Eq. 
These in turn determine the allowed phases, the ordering temperatures and associated wavevectors, and the relative 
magnitudes of the relevant order parameters in mixed states. We considered several other crossings in less detail than 
the one presented here, and these investigations, combined with general expectations, allow us to outline the features 
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FIG. 6: Phase diagram at vt 



8.0 in the q — T plane. Temperature is in units of 




and a is dimensionless. The 



degeneracy point is indicated by a dotted black line on the a axis. The shaded region indicates where SSE ordering is preferred 
over TTE ordering. Temperature Tzz is below the temperature range. 

we believe to be more or less generic for CDW ordering when there are two levels close to degeneracy. The vt-T phase 
diagram we show in Fig.|3]is qualitatively very similar to the temperature-filling phase diagrams found in the case of a 
single LL by Fogler et al. (FKS)i^ and Moessner and Chalker (MC)i^ if we divide = vi + v-i by 2. In both of these 
works it was found that for v close to \ in the valence LL, there is a striped CDW, with a transition to a triangular 
CDW at zero temperature for > 0.35 (FKS) or > 0.39 (MC). Taking a similar approach, we extrapolate the 
line separating the SSE and TTE phases in Fig. U to zero temperature and find that the SSE phase is stable when 
the effective filling fraction of each state is approximately fc > 0.45. As MC point out, the perturbative expansion of 
the free energy is not valid at very low temperatures, so our estimate of Vc is only approximate. Additionally, we can 
reasonably expect that the precise estimate of i>c may depend on which two levels are being studied. In this sense, 
the results here for degenerate energy levels are consistent with known results on single-LL systems. However, there 
is an important difference, in that there are effectively two half-filled states and there is co-ordinated ordering in both 
of them, as opposed to the single half- filled Landau level studied by FKS and MC. 

Away from degeneracy the TTE phase tends to be dominant. As vq increases towards unity, the SSE phase has a 
higher transition temperature just near degeneracy. In all of the phases where there is ordering in both states 1 and 
2, we observed that the wavevectors Qi = Q2, which can be understood as arising from the extra free energy that 
can be gained when the wavevectors are equal in the two states and the order parameters are in phase, as illustrated 
in Eq. 

At degeneracy, we compared the transition temperatures of the CDW states to the QHF states, and found that the 
system preferred CDW ordering over the QHF ordering. We are uncertain as to the extent to which this behaviour 
is generic, however, our results appear to indicate that the CDW state can be competitive with quantum Hall 
ferromagnetism. 



We have considered charge density wave ordering in a 2D hole system in a perpendicular magnetic field when 
Landau level mixing is important. The motivation for this work is recent transport experiments in quantum Hall 
systems where Landau level mixing was induced by spin-orbit coupling i^^i^^ Our work is related to the situations 
investigated by Manfra et al~- and Fischer et alM There has been some previous theoretical work on both hole 
striped states^ and hole FQHE states,^® but these works do not focus on the situation when two levels are close to 
degenerate and the possible orderings in such a case, as we do here. 

Recently, Manfra et al^ studied the transport properties of a two-dimensional hole system in a perpendicular 
magnetic field. They found anisotropic resistivity at some, but not all, half-integer filling fractions; in particular, they 
observed isotropic charge transport at 1^—9/2, flanked by anisotropic transport at iy—7/2 and i/=ll/2. They argue that 
LL mixing caused by spin-orbit coupling in the valence band of GaAs may be responsible for this pattern of isotropic 
and anisotropic transport. From a self-consistent calculation of the Landau levels in the Hartree approximation, they 
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suggested that the contribution from the = and = 1 LLs to the wavefunction of the v — 9/2 state in their 
particular device hkely suppresses the tendency towards charge density wave order. 

In Fig. [5] we study the situation in which j/q = 0.5 and find that weU away from degeneracy so long as there is some 
occupation of the upper energy level i^i ^ 7^ the effect of LL mixing can be to stabilize a mixed triangular CDW 
state over the striped CDW state we would expect for an unmixed Landau level. This is another mechanism that 
could lead to the absence of anisotropic transport at = 0.5. 

Fischer et al^ observed anomalous transport in the lowest Landau level at what they interpreted to be a field- 
induced level anti-crossing. At a temperature of 320 mK there was a peak in the longitudinal resistivity that disap- 
peared at 50 mK. As we do not expect our theory to be applicable to the lowest Landau level, we do not speculate 
on the origin of this anomalous transport. However, we do note that their experiment exhibits a level anti-crossing 
that occurs because they study a two-dimensional hole rather than a two dimensional electron system. Rashba and 
Dresselhaus coupling play an important role in determining the energy levels in their device, and if they tuned the 
Rashba coupling, they would be able to explore the effects of tuning mixing at = 0.5 and = 1 directly. 

Our calculations generalize existing work on 2DEGs in high-LLs to the situation when LL mixing is important, 
and we show that for a contact potential the Hartree-Fock approximation is exact in the high LL limit even in the 
presence of LL mixing. We specialize to the situation in which two levels are near the Fermi energy and perform a 
diagrammatic analysis of the two-particle vertex. This leads us to predict three relevant temperature scales where 
the uniform system becomes unstable to charge-density wave formation. Two of these correspond to CDW formation 
in either of the two levels, and the third corresponds to an instability to a mixed CDW phase. Having established a 
tendency towards CDW ordering, we derive a Landau theory for CDW ordering in the two levels using the Hartree- 
Fock approximation. We then focus on the specific example of LL mixing induced by Rashba spin-orbit coupling 
(which could provide a tuning parameter for experimental investigations of the effects of LL mixing on anisotropic 
transport in the quantum Hall regime). 

The effect of Landau level mixing goes beyond changing the character of the single-particle states in the system to 
also changing the spectrum and the spacing in energy of the single-particle states. We find that it is the second of 
these that appears to have the most effect on anisotropic quantum Hall states, once the single-particle states have 
little A^ = or 1 character. In particular, at the level of Hartree-Fock, if there are two states relatively near the Fermi 
energy, only when the filling of one or both states is very close to v = 0.5 is there striped CDW ordering. At other 
fillings there is triangular CDW ordering. A competing phase when there are two close to degenerate states is that 
of quantum Hall ferromagnetism, and while in the example we consider here we find the ordering temperature to be 
less than that for a CDW, thermal fluctuations may treat the two states differently. 

We hope that our work stimulates further theoretical and experimental study of anisotropic Quantum Hall states 
in which externally tuned LL mixing is used as a parameter to investigate the phase diagram. 
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In this Appendix, we give exact expressions for the coefficients in the Landau free energy theory for each case 
discussed in Sec. IIIII when Qi 7^ Q2. In solving for the Ci and 7 terms, we determine the chemical potential by setting 
the number of fermions in the non-interacting system equal to the number of fermions in the interacting system. 
Then, 
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Appendix A: Full expressions for coefBcients in Landau theory 



1 Tkk' E{Q,} n Wkk']" (1 - ^fcO (1 - 2^fe0 



(Al) 



1. Mixed unidirectional-unidirectional CDW system 



The general form of the free energy for striped phases is 



:F = Pfai (Qi, T) |Ai|' + P|a2 (Q2,r) IA2I' + F^i (Qi,r) |Ai|* + P^c^ (Q2,T) IA2I' 
+F2p2^(Qi,Q2,^,r)|AinA2^ 



(A2) 
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Note that terms involving A? are not allowed by symmetry. The parameters for the free energy theory in the case of 
mixed unidirectional-unidirectional CDW ordering are 



ais {Qi, T) = Uii {Qi) + Uu {Qif f {Ei - m) + Uij {Qif f {Ej - m) , 



(A3) 



Cis{Qi,T) = 



and 



7ss(Qi,Q2,r) 



Uii {Qif f" [Ei - + ^ Uij {Qif r {Ej - fj,) 



; { Uu {Qif f" {Ei - m)' + E [^ii (^')' - 



+ 2Uii {Qif Uij {Qj) f" {Ei - /i) /" {Ej - m)] I , 



{Qi)Uu {Q2ff" {Ei-ii) 



+ \Ul2 {QlfU22 {Q2f f"'{E2-lj) 



2 + cos {Q1Q2 sin 6*) 
2 + cos(0iQ2 sin (9) 



1 



/" (El - m) /" {E2 - m) Uii {Qif U22 {Q2f 



Y.if'{Ei-n) 

+f" (El - m)' U^i {Qif U12 {Q2f + f" {E2 - m)' Ui2 {Qif U22 {Q2f 
+f" (El - /" (E2 - /x) C/12 (Qi)' U21 {Q2f^ 



(A4) 



(A5) 



2. Mixed triangular-triangular CDW 

The free energy when there is mixed triangular-triangular ordering is 

= P^ai (Qi,T) I All' + P|a2 (Q2,T) IA2I' + Pfh (Qi,T) |Ai|' + P|62 (Qi,r) IA2I 
+ PiV (Qi,T) |Ai|* + Pic2 (Q2,T) IA2I' + PfP|7(Qi,Q2,^,T) I All' IA2I' , 



(A6) 



A triangular CDW is described by a set of three wavevectors | Q^, | lying on a circle of radius Qk oriented ^ radians 

apart. Allowing for non-coUinear CDW ordering between two triangular CDW states, we define 6 as the angle between 

wavevectors Qj and Q2. Then, the angle between any two wavevectors in $1 and $2 is ^Qi, Q2 ^ = ^ + ^ (i' ~ j) ■ 

There are three possible angles between wavevectors: {6i}^^i = {9,6 + 6 — 
The parameters in this free energy model are 



aiT {Qi,T) 



U^i {Qi) + Uu {Qif f (E, - /x) + E C/y {Qif f (E, - 11) 



(A7) 



h{Q^,T) 



-2 COS 



Uu {Q^f f" (E. - + E (Q^)' /" (^^- - ^) 



(A8) 
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1 — cos 



9 



2E. /' (^.-m) 



U^^ {Q^t f" {E, - flf 



E 



/" - (g,) + 2/" - /" {Ej - ii) t/,, (q,) L/y (q.) 



and 



7TT (Qi, Q2, T) = f/n (Qi)' [/i2 (Q2)' /"' (-El - /i) + (Qi)' C/12 (Q2)' /"' (-^i - /^) 



X 6 + cos {Q1Q2 sin 0) + cos 
9 



2tt 

]iQ2Sm[e + — 



cos 



, 2tt 
Q1Q2 sm{9- — 



/" - /i) /" (i?2 - /i) f/ll (Ql) t^22 (Q2)' 



+/" - Ai)' Un {Qif U12 (Q2)' + /" iE2 - fif U12 (Qif U22 (Q2)' 
+/" - fi) J" (E2 ~ fi) C/12 {Qif U22 (Q2)' 



(A9) 



(AID) 



3. Mixed unidirectional-triangular CDW 

Without loss of generality, assume that there is unidirectional CDW ordering in $1 and triangular CDW ordering 
in $2- Then the free energy is 



F = p2^i(Qi,T)|Ai|VP|a2(Q2,r)|A2|' + P|fe2(Qi,T)|A2|' 

+P4ci (Qi,r) |Ai|Vp|c2 (Q2,T) IA2I* + P2p| ^(Q,,Q2,0,T) jAij^ IA2I 

All coefficients except 7 remain unchanged for each type of ordering. The coupling parameter is 
1ST (Qi, Q2, T) = i [c/n {Qif C/12 (Q2)' /"' {El - m) + C/n (Qi)' C/12 (Q2)' /"' - m) 



X 6 + cos {Q1Q2 sin 6') + cos 
3 



2tt 

)i(92sin [0+Y 



cos 



Qi(52 sin — 



/" {El - f" {E2 - /i) c/11 {QiY U22 {Q2r 



+f" {El - ^lf Uii {Qif c/12 (Q2)' + /" {E2 - fJ.f U12 {Qif U22 (Q2)' 
+/" {El - /i) /" {E2 - c/12 (Qi)' U22 (O2)' 
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Appendix B: Full expressions for coefficients in Landau theory when Qi = Q2 



In the previous appendix we gave exact expressions for the coefficients in the Landau free energy theory for each 
case discussed in Sec. IIIIl When the ordering wavevectors Qi = Q2 = Q then there are additional terms in the free 
energy, as described in Eq. as <5-7^Qi=Q2- We give expressions for the coefficients in Eq. ([M)) in this Appendix. 
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1. Mixed unidirectional- unidirectional CDW system 



The extra terms that appear in the free energy theory for SSE ordering are of the form 

SfssE = PiP2assE{Q,S(f>,T)\A,\\A2\ + P^P2assE{Q,5<f>,T)\A,f \A2\ 

+PiPhssE {Q,5ci>,T) |Ai|2 IA2I' + PiPipssE {Q,5ci>,T) |Ai| IA2I' . 



(Bl) 



Hero, Q is the magnitude of the ordering wavcvcctor, 54) the relative phase between order parameters Ai (Q) and 
A2 (Q), and T the temperature. The values of the coefficients are 



assE{Q,5(l>,T) = 2Ui2{Q) 



1 + f/ii (Q) /' {El - + U22 {Q) f {E2 - lA 



cos {5(j)) , 



(B2) 



CTSSE (g, Sep, T) = { Un {Qf t/12 (0) /"' {El - fi) + U12 {Qf (Q) /"' {E2 - fi) 



, Uii{QY f" {El- n) + Ui2{QY f" {E2- IJi) 
X [Un {Q) U12 (Q) /" {El -n) + U22 (g) t/12 (g) /" {E2 - fi)] ] COS , (B3) 



7ssE{Q,S(t>,T) = <i- 



Un {Qf U12 {Qf f" {El -n) + U12 {Qf U22 {Qf f" {E2 - IJ) 



T.Uf'{Ei-^^) 



X [Un (g) U12 [Q) /" {El - /x) + U22 (g) U12 (g) /" {E2 -^l)f\[l + 2 cos" {5(t>)\ 



Eti /' {Ei - n) 



Un {Q) U22 {Q) - U12 {Qf f" {El - fi) f" {E2 - m) , 



and 



(B4) 



PssE{Q,Sct>,T) = {Un{Q)Ui2{Qff"{Ei-f,) + Ui2{Q)U22{QTf"'{E2-tJ-) 



3 fi/i . 



U22 {Qf f" {E2 - m) + U12 {Qf f" {El - n) 



T.Uif'{Ei-i^) 

X [Un {Q) U12 {Q) f" {El - m) + U22 {Q) U12 {Q) f" {E2 - cos {dcj>) . (B5) 



2. Mixed triangular-triangular CDW system 

The order parameter for the triangular CDW state depends on a set of three wavcvcctors, {Qf}^^i- We associate 
ordering at each wavevector with a particular phase (pf, such that Ai (Qf ) — A^e*'^". Thus, when Qi = Q2, we end 
up with six new parameters to consider when minimizing the free energy. We may reduce the number of parameters 
from six to four in the following manner (where i = 1, 2, 3): 



Vl - V2^ 
i=l 



(B6) 
(B7) 

(B8) 
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The new terms in the free energy theory for TTE ordering are 

5fTTE = PlP2aTTE {Q,Sct)l,5(^2,Sct>3,T)\Ai\ \A2\+ P^P2VTTE (Q, 5-^1, 5-^2, <5</'3, ^1, t) |Ai|^|A2| 

+PlP^ipTTE [Q, 5cj)i, Sct>2, Sct>3, 4>2,T) |Ai| |A2|^ + P^P2<JTTE (Q, ^<^1, 5<^2, Scj)^, T) |Ai|^ I A2I 

+PlPhTTE iQ,ScPi,5cj>2,6cj>3,T) I All' IA2I' + PiPipTTE {Q,5(^i,5(t>2,5h.T) I All IA2I' . (B9) 

Here, Q is the magnitude of the ordering wavevector, dcj) the relative phase between order parameters Ai (Q) and 
A2 (Q), and T the temperature. The values of the coefficients are 



axTE {Q,S(/)i, 6(1)2, S(/)3,T) = Ui2{Q) 



1 + Uii (Q) f (£1 - /z) + U22 (Q) f {E2 - /x) 



E 



cos 



(BIO) 



Vtte (Q,S(t)i,5(j)2,5(j)3,4>i,T^ = 



-6 



f/ii {QY f/i2 (Q) /" (El - m) + Un (Q) f/i2 {QY f" {E2 - 



X COS I j ^ cos - 6(j)i^ 



(Bll) 



Vtte (Q,S(j)i, 6(1)2, 5(p3,4>i,T^ = 



-6 



U12 {qY Uii (Q) /" (£1 - m) + U12 (g) U22 {qY f" {E2 - lA 



X cos I j ^ cos (^02 + , 



(B12) 



i=l 



ctttb {Q,6(I),T) 



11 + 4cos(^Q2 



{QY U12 iQ) f" {El - m) + U22 (Q) Ui2 {QY f" {E2 - m) 



E-=i/'(^^i-M) 



[Uii iQ) C/12 (Q) /" {El - ^) + C/22 (Q) [/12 (Q) /" (J^2 - m)]^ 



Un {QY + U12 {QY 



COs((?!>ia -(j)2a): 



(B13) 



7TTB {Q,6(j),T) 



Un {QY U12 {QY /"' (£^1 - + U12 {QY U22 {QY f" {E2 - n) 
1 



Eti /' {E^ - m) 

3 

x^[l + 2cos2((5(/)i)] 

i=l 

'2 " 



[f/ii (Q) U12 {Q) f" {El - /x) + U22 (Q) «7i2 (Q) /" {E2 - 



c/11 (g)' U12 {Qf /"' (^1 - m) + U12 {Qf C/22 (g)' /"' {E2 - fi) 



2 + cos(^Q'^ 



Eti /' {Ei - m) 
X [Un {Q) U12 {Q) f" {El - + U22 {Q) U12 {Q) f" (£2 - 



X ^ ^ [1 + 2 cos {6(j)i) cos 
9 



(g) U22 (g) - c/12 (g)' /" {Ei - n) f {e^ - n) , 



(B14) 



23 



and 



Ptte {Q,H,T) 



11 + 4 cos i^Q^ 



Ui2 {Qf U22 (Q) /"' {E2 + C/11 (Q) U12 {Qf f" {El - 11) 
[f/11 (Q) {Q) f" {El - + U22 (Q) C/12 {Q) f" {E2 - n)] 



U12 {Qf + U22 {Qf 



^ cos {(j) 



la — (P2a) 



(B15) 



3. Mixed unidirectional-triangular CDW system 

We consider mixing between a unidirectional phase in $1 and a triangular phase in $2- Without loss of generality, 
we let Qi = Q2- Thus, we end up with two new parameters: 

Sct> = (t>i-(l>2, (B16) 
3 

^2 = (B17) 

i=l 

The new terms in the free energy theory for STE ordering are 

SfsTE = PiP2asTE {Q,5(I>,T) I All IA2I + PiPi^sTE {q,5(I>,^2,t) I All IA2I' + P^PhsTE {Q,S(t>,T) I All' IA2I' 
+PiP2(TSTE {Q,S4>,T) \A1fA2 + PlPipTTE {Q,S4>,T) |Ai| |A2|^ . (B18) 
The coefficients are 



asTE{Q,5cj>,T) = 2Ui2{Q) 



1 + f/ii (Q) /' {El - /x) + U22 {Q) f {E2 - m) 



cos {S(j)) . 



(B19) 



<fiSTE(Q,S<P,^2,T'^ = -2 



U12 {Qf Uii (Q) /" {El - + U12 {Q) U22 {Qf f" {E2 - 



X cos ^"^Q^^ COS {4>2 + S(jfj , 



(B20) 



o-STB {Q,S(I),T) 



19 + 8 COS (^Q2 



Uii {Qf U12 (Q) /"' (El + U22 {Q) C/12 {Qf /"' {E2 - m) 
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E^=l f {E^ - Ai) 

X (Uii {Qf + U12 {Qf] ] cos {S<P) , 



(B21) 



ISTE {Q, S^, T) = {^Uii {Qf U12 {Qf f" {El - m) + U12 {Qf U22 {Qf f" {E2 - fJ,) 

[Un (Q) C/12 (Q) /" {El -n) + U22 (Q) C/12 (0) /" {E2 - fi)f 



Y:Uf'{Ei-f^) 



+ 



1 + - cos' {6(1)) 



3 + 2cos[^Q' 



E-=i/'(^^-m) 
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[Un (Q) C/12 (0) /" {El - m) + C/22 (Q) C/12 (0) /" (i52 - p)f \, (B22) 
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and 



PSTE {Q,6(t>,T) = I - 



19 + 8 COS [^Q2 



Ui2 (Qf U22 (Q) f" {E2 - + [Q] f/i2 {Qf f" {El - m) 
[Uii iQ) Ui2 (Q) /" {El + U22 (0) Ui2 {Q) f" {E2 - fi)] 



Ui2 {QY + U22 {QY 
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Appendix C: Hartree-Fock potentials for quantum Hall ferromagnets 

In this appendix, we derive expressions for the potentials Wij in the Hartree-Fock approximation. We use notation 
similar to Sec. IIIIl but note that lower-case letters label the individual states while upper-case letters label the 
pseudospin indices. Let 



(9^0)' 



where 



Then 



p(q) 



'-x"^ C-'^ {x) , m<n. 



We invert Eq. ([57)) to solve for o,jX- in terms of the rhj^q operators and rewrite p (q) as follows: 



P(q) = e 4 ^ J"/(q)TO/, 

IJ—0,x,y,z 



where J-j (q) — ^a'"/'^^ Fij (q). The density-density Hamiltonian is then 



(CI) 



Hqhf = ^^^ij{^)^i,cim,j-ci, 



q IJ 



(C2) 



where Vjj (q) — v {q)J-i {<l)J-j (— q), v (q) being the electron-electron interaction potential. 

We apply the HF approximation to Eq. (IC2p to obtain Eq. (15^]) . The HE potentials Wjj are defined as follows: 



Woo = 


2^ j 


f d\e 


Wii = 


2^ j 


fd\e 


Wij = 


2^ j 


f d\e 



Voo (q)-^oo (O)--u(q) 
Vii {q)~Vii {0) + lu{q) 



where 



u(q) = Fu(q)F22(q)-^^i2(q)i^2i(q). 



Eor each potential, conservation of LL index in the basis state interactions must be observed. Thus, for the cigcnstates 
considered in Sec. |lVl the only nonzero potentials are W^oo, Wzz, Woz — W^Oi and Wxx — Wyy. 
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